Subgap noise of a superconductor-normal-metal tunnel interface 
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It is well established that the subgap conductivity through a normal-metal-insulator- 
superconductor (NIS) tunnel junction is strongly affected by interference of electron waves scattered 
by impurities. In this paper we investigate how the same phenomenon affects the low frequency 
current noise, S, for voltages V and temperatures T much smaller than the superconducting gap. If 
the normal metal is at equilibrium we find that the simple relation S(V,T) = 4ecoih(eV/T)I(V,T) 
holds quite generally even for non-linear I-V characteristics. Only when the normal metal is out of 
equilibrium, noise and current become independent. Their ratio, the Fano factor, depends then on 
the details of the layout. 
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PACS numbers: 74.40.+k, 74.45. +c, 74.50.+r 



I. INTRODUCTION 

Recently, a great progress has been achieved in 
the theoretical understanding of current fluctuations in 

Smesoscopic normal-metal superconductor (N-S) hybrid 
systemsi^ This progress has been boosted by the devel- 
opment of simple techniques to calculate the full count- 
ing statistics of quantum charge transfer. 4-5-6 In partic- 
ular, the current fluctuations in a diffusive wire in good 
contact with a superconductor have been calculated tak- 
ing into account the proximity effect for any voltage and 
temperature below the gap^ With proximity effect we 
mean here the presence of a space dependent coherent 
propagation in the normal metal of electrons originating 
from the superconductor. The opposite limit of a tunnel 
junction between a diffusive metal and a superconductor 
has been less investigated, partially because only very re- 
cently it has become possible to measure current noise in 
tunnel junctions^ Theoretically, current noise in a NIS 
junction has been considered by Khlus long ago^ but ne- 
glecting proximity effect. Later, de Jong and Beenakker 
included the proximity effect at vanishing temperature 
and voltage^ The effect of a finite voltage was studied 
very recently using a numerical approach. 11 More com- 
plicated structures with several tunnel barriers have also 
been considered. In some limits these can be reduced to 
a single dominating NIS junction with a complex normal 
region^2ii4 Actually, one may expect that the noise at 
finite voltage and temperature (of the order of the Thou- 
less energy E t h) depends on the spatial layout when a 
tunneling barrier is present^* but this has not been in- 
vestigated so far on general grounds. 

As a matter of fact, it has been shown by one of the au- 
thors and Nazarov in Ref. 

El 

that the subgap Andreev 
tunnel current is strongly affected by the coherent scat- 
tering of electrons by impurities near the junction region. 
Two electrons originating from the superconductor with 
a difference in energy of e can propagate on a length scale 
of the order of £ e = y/D/e before dephasing (D being the 
diffusion coefficient and we set K and ks = 1 throughout 



the paper). If the relevant energy scale of the problem, 
i.e. the voltage bias V multiplied by the electron charge e 
and the temperature T are sufficiently small, the result- 
ing coherence length £ corr = \J Z?/max{eV, T} is much 
larger than the mean free path I. Thus at low temper- 
atures and voltages the electron pairs are able to "see" 
the spatial layout on a length scale given by £ corr ^> I- 
Since electron pairs attempt many times to jump into the 
superconductor before leaving the junction region, inter- 
ference enormously increases the current at low voltage 
bias and the resulting conductivity depends strongly on 
the explicit layout^ In this paper we investigate how the 
same coherent diffusion of pairs of electrons determines 
the subgap current noise of a NIS tunnel junction. 

We find that a generalized Schottky relation holds for 
voltages and temperatures smaller than the supercon- 
ducting gap: 



S(V, T) = 4e coth(eV/T)7(V, T). 



(1) 



Due to the interplay between the proximity effect and the 
presence of the barrier the current voltage characteristics 
is both non-universal and non-linear. However, accord- 
ing to Eq. 0, the ratio F = S/(2eI) = 2coth(eV/T), 
known as Fano factor, is universal as long as the elec- 
tron distribution on the normal metal is at equilibrium. 
In particular, for T <C eV the Fano factor is 2, indicat- 
ing that the elementary charge transfer is achieved by 
Cooper pairs. 

If the normal metal is driven out of equilibrium, within 
our formalism F can be calculated once the geometry is 
known. One finds then that the noise and the current 
remain both independent of each other and strongly de- 
pendent on the geometry. We discuss a realistic example 
(cf. Ref. [H where we predict strong deviations from 
Eq. HJ. 

The paper is organized as follows. In Section [H] we 
set up the main equations and obtain the result given 
in Eq. 0. In Section ITTTI we discuss the case when the 
normal metallic reservoirs are out of equilibrium. In this 
case Eq. does not hold and to obtain the Fano factor 
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the diffusive nature of transport will play an important 
role. Section HVI gives our conclusions. 



II. TUNNELLING HAMILTONIAN APPROACH 

Let us consider a normal metallic reservoir connected 
through a tunnel junction to a superconductor. This 
problem can be treated conveniently by perturbation the- 
ory in the tunnelling amplitude with the standard model 
Hamiltonian H = Hs + + Ht- Here Hm and Hs de- 
scribe the disordered normal metal and superconductor: 



Hn = £64 CT e kCT , (2) 
Hs = EC k <V + E[ A <l d -qT+ A * d - q T d q lK?) 



A is the superconducting order paramenter, Cko- and d qcr 
are destruction operators for the electrons on the nor- 
mal and superconducting side, respectively. With k and 
q we indicate the eigenstates of the disordered Hamilto- 
nian with eigenvalues £k and £ q , q and — q indicate time 
reversed states. The index a stands for the spin projec- 
tion eigenvalue. The tunnelling part of the Hamiltonian 
reads: 



ht=J2 fa 

k,q,cr 



(4) 



where ikq = / drdr't(r, r')-0k(r)'0q( r/ )j and t(r, r') is the 
quantum amplitude for tunnelling from point r of the 
normal side to point r' of the superconducting one. All 
the information on the disorder is in the eigenfunctions 
V'k/q- The difference of potential V between the two elec- 
trodes can be taken into account by the standard gauge 
transformation £k £k — e.V . 

We are concerned only with sub-gap properties: 
eV, T <C A. No single particle states are available in 
the superconductor for times longer than 1/A. Thus the 
low frequency noise (to <C A) will be determined only 
by tunneling of pairs. Single particle thermal excitations 
in the superconductor are exponentially suppressed and 
are thus neglected. To calculate the quantum amplitude 
for transferring one Cooper pair to the normal metal we 
proceed as in Ref. 0. This gives: 



A k k' = ^t kp t k /_ p u p fp 



, (5) 



where k and k' indicate the final electron states (assumed 
empty), up = 1 - u\ = (1 - C P /£ p )/2 are the BCS co- 
herence factors, and E p = (£ 2 + A 2 ) 1 / 2 is the supercon- 
ducting spectrum. 

The function A gives the amplitude for the only pos- 
sible elementary process at low energies. It is thus con- 
venient to write an effective tunnelling Hamiltonian in 



terms of these processes, 
H eff =H N + J + J^ with J = J2 A ^' c U c li b o, 

kk' 

(6) 

where b = (2/Ns) X) q 0q^qT^-qi i s the destruction op- 
erator for a Cooper pair in the condensate, q = w q /w q is 
the pair wavefunction and N$ is the number of fermions 
on the superconducting side. Note that b is not a true 
bosonic operator, since Cooper pairs overlap— In prac- 
tice, since the superconductor is in the coherent phase, 
the averages involving b are trivial to perform (we ne- 
glect collective-modes excitations, which are at higher en- 
ergy): (6 ) = {b\) = (b\b ) = I. A similar Hamiltonian, 
but with a constant amplitude, has been considered by 
many authors to study the influence of collective modes 
on Andreev reflection from a phenomenological point of 
viewiiSii 9 . Here we found its form from the BCS micro- 
scopic model, the resulting H e H is thus equivalent to the 
starting one (0J at low energies. 

The current operator— can be obtained as usual by 
the time derivative of N = X) qCT c Lr c kcr : 



dN 



Kt) = e =1 = ei[H°", N] = -2ei[J(t) - j\t)} . (7) 

We have now all the elements to calculate by perturbative 
expansion in the new coupling A (~ t 2 ) the quantum av- 
erage of the current, (/), and of the current-noise-spectral 
density: 



S(u = 0) 



dt[([I(t) ) I(0)] + )-2(I(t))(I(0))]. (8) 



Writing the evolution operator in the interaction rep- 
resentation U(t) = Tcxp{-i f ^dt'iJoit') + 4 (*')]}, 
where J (t) evolves according to H^, one can readily 
expand in A and obtain in lowest order the following ex- 
pression for I and S: 



I{t = Q) = 2e[N^-N^} 
S(u = 0) = 8e 2 [N^ + N^' 



(9) 
(10) 



where 



N^= dt'(J o (0)4(t')) , N^= dt'(4(t')Jo(0)) 



This result depends only on the fact that the current 
can be described as a first order process in the effective 
tunnelling amplitude A. For the tunnelling case Levitov 
and ReznikoAiSS showed that the full counting statistics 
is bi-poissonian: 



lnx(A) = (e 2leA - 1)JV«_ + (e~ 2leA - 1)JV. 



-2ie\ 



(11) 



where x is the generating function. Taking successive 
derivatives with respect to iX of In \ and setting A = 
afterwards one generates current momenta of all or- 
ders. It is readily verified that / = <91nx/<9(iA), and 



S = 2d 2 \n X /d(iX) 2 coincide with ©-Jig. Eq. fH} im- 
plies that / and S determine the full counting statistics 
of the problem. The physical picture is simple: the dis- 
tribution of transmitted Cooper pairs is given by a su- 
perposition of two Poissonians for tunnelling from left to 
right and from right to left. 

If the metallic side is at the thermodynamic equilib- 
rium the full dependence of noise on voltage and temper- 
ature is given on general grounds by ft}. This relation 
can be proven by evaluating © and (fTUf) using the basis 
of eigenstates \n) of with eigenvalues E n . One can 
actually write / = 2eC_ and S = 8e 2 C+ with 



C± = 



n m 



e 



dt „2ieVt 



,-E n /T 



\2 iuj nrn t 



i I J nm I ^ ) 



(12) 



where uj nm = E n — E m , Z = J2 n e~ En ^ T is the partition 
function and J nm — (n\J\m). By simple manipulations 
of Eq. 1|12[) one obtains 



C± 



2tt( 



n rn 



2 5{2eV - uo nm ) . 

(13) 

From Eq. (|13|l the relation JIJ follows directly. Note that 
the proof is valid for any system at thermal equilibrium, 
interactions in the normal metal do not spoil the result. 
For a generic tunnelling system the relation between noise 
and current was shown by Sukhorukov and Loss^i and 
discussed recently by LevitoviSS 

One thus finds that non-linear dependence on the volt- 
age and temperature of the current will be found exactly 
in the low frequency noise with a universal Fano factor 
of F = S/(Ie) = 2coth(e^/T). In the completely coher- 
ent case (^4 independent of energy, i.e. eV <C E t h) the 
same equation has been proved in Ref. |U In this paper 
we extend its validity to arbitrary ratio eV/ E t h provided 
eV <C A. Few comments are in order. Eq. can be re- 
garded as a generalization of the fluctuation dissipation 
theorem, since it holds for any value of the perturbing 
field V~ We would point out also that in presence of 
interactions among electrons the noise was calculated re- 
cently in Ref. with the limitation eV <C E t h- We men- 
tion that for a normal wire in good contact with the su- 
perconductor the situation is much richer and the full en- 
ergy dependence of the Fano factor has been calculated 7 
and measurediSiiSSiS 7 . For this case an analytical expres- 
sion that links the current to the noise was proposed, 28 
but it does not reproduce the weak energy dependence of 
the differential Fano factor [cf. Eq. (|28|l in the following] . 
Only very recently an analytical solution of the gener- 
alized Usadel equation introduced in Ref. has been 
obtained^ 

In order to emphasize the generalityof Eq. let us 
consider a non-trivial example. In Ref.|3fjthe low- voltage 
anomalies due to the proximity effect have been de- 
tected in a superconducting ring closed on a normal metal 
through two tunnel junctions as shown in Fig. (Q. Part 
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FIG. 1: Experimental set up similar to Ref l30t The mag- 
netic flux dependence of the current and noise can be used to 
measure the doubling of the Fano factor with high accuracy. 



of the subgap current through this structure is modulated 
by the external magnetic field. Using Eq. JQ) one can pre- 
dict the following noise dependence on the magnetic flux 
$ through the ring: S($) = 4e coth(eI//T)/($). Since 
the <E> dependent part is given by the coherent propaga- 
tion of two electrons between the two junctions any single 
particle contribution is eliminated from the outset. This 
experiment could thus be used to test accurately the va- 
lidity of Eq. Q and in particular the doubling of the 
charge at low temperature. 



III. OUT OF EQUILIBRIUM RESERVOIRS IN 
DIFFUSIVE JUNCTIONS 

When the normal metal is not at equilibrium the sim- 
ple relation Q between S and / does not hold any longer, 
and expressions © and IjlUI) have to be evaluated explic- 
itly. We proceed by writing: 



/ = 47re^|A kk f + 2eV)Jf_(&,&0 , 

kk' 

S = 167re 2 ^|A kk f 6{& + & - 2eV)fl+ 



kk' 



where: 

H±& O = [1 - n(0] [1 - n(€')] ± n(0»(O . (W) 

and n(£k) = ( c l ia c ka)- ^ ls convenient to introduce the 
relevant energy dependence by defining: 

A(e,eV) =^5{e/2 + eV-&)5{-e/2 + eV-&)\Aw\ 2 . 

kk' 

We thus obtain: 



, = 2vre/ deAle, eV) . 



(15) 
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If the normal metal is in thermal equilibrium is the 
Fermi distribution f F (Q- In that case H + and if_ are 
related by the simple relation: 



H+(e/2 + eV,-e/2 + eV) , ( eV 

; ; ; = COth 

H_(e/2 + eV,-e/2 + eV) \T 



(16) 



and we recover again Eq. QJ. If ^ /f(£), Eq. I|16|) 
does not hold and S and / become independent. Indeed, 
the Fano factor ceases to be universal and gives new in- 
formation on the system. The technique developed so far 
enables calculation of F also in this case through Eq. (|15fl . 
however now we need to specify the explicit form of A(e). 
We follow Refs. to evaluate it: 

A(e,eV) = J dC<«*XC;£,r:W;£,0>< 
S(C, C; C, e)S(e/2 + eV- £)8(e/2 - eV + £') 

where Ffot,?) = u(CMC){(£ + - £(0)^ + (f + 
eV-EiOr 1 }, 

S[C,C',e,e'] = yd 3 ri...rf 3 r 4 Jd 3 r' 1 ...d 3 r' 4 x 



t*(r lj ri)t*(r 2 y 2 )t(r3 ) ^)t(r 4j r / 4 ). 



(17) 



Merely tf c (ri,r 2 ) = Ek ^(C - Wk(ri)^(r 2 ) is the 
single particle propagator and (. . )i mv indicates impu- 
rity average. Since the tunnelling is most probable only 
near the surface of the junctions, we can set = and 
restrict them to lie on the surface. The impurity aver- 
ages can be performed within the usual approximation 
kpl 3> 1 (k F Fermi momentum). The main contribu- 
tion to Eq. 117fl is given by the "Cooperon" diagrams, 
associated to the coherent propagation of two interfering 
electrons with small energy difference e and small total 
momentum. The range of the Cooperon propagator is 
£ e , and thus it is much larger than the range of K^(r, r') 
which is I. For this reason among the integrals over cPr 
two dominant terms can be singled out: (a) the term 
where ri as r 2 and r3 as r 4 with interference occurring in 
the normal metal and (b) the term where r 4 as r3 and 
r 2 as r 4 with interference occurring in the superconduc- 
tor. Since in Ref. it has been shown that the second 
term is always smaller than the first by a factor eV/A we 
consider only case (a). The dominant contribution to S 
can thus be written as follows: £, £'] as S(£ — £') 

and 



S(e) 



32n 3 e 4 v N S 2 



J rf 2 r 1 d 2 r 2 [P e (r 1 ,r 2 )+F_ £ (ri,r 2 ) 



_(18) 

Here Gt is the conductance in the normal phase, S is the 
surface area of the junction, vn is density of states of the 
normal metal (per spin) and P e (r) satisfies the diffusion 
equation in the normal metal: 



The integration over £ and £' for £, £' ,eV <C A can be 
done and gives 7r 2 , we are thus left with an explicit ex- 
pression for A(e) = 7r 2 S(e) which can be calculated if the 
geometry of the system is known. Given A and the 
noise and the current can be found explicitly from H15|) 
and UIJ. 



A. An explicit example: a wire out of equilibrium 

We are now in a position to calculate both current 
and noise in a non-equilibrium system. Let us consider a 
realistic example: In Ref.ll6lthe non-equilibrium electron 
distribution for a small metallic wire of length L has been 
measured by using a tunnel junction between the wire (at 
different positions x along the wire) and a supercondutor 
(see also inset of Fig. [3] in the following). In that case 
the quasiparticle current was measured, but the subgap 
current and noise in the same configuration could also 
be measured. When L is shorter than the inelastic mean 
free path the electron distribution along the wire is given 
by the diffusion equation*^ 

n(x, = f F (£ + eU)[l - x/L] + f F {£)x/L (20) 

where U is the difference of potential between the two 
normal reservoirs. This leads to a double discontinuity 
of the distribution function at T = 0. 

Since the transverse dimension w of the wire is much 
smaller than £ CO rr, we can use the one dimensional diffu- 
sion equation. Following Ref. |^ we impose the bound- 
ary condition on the propagator P corresponding to a 
good contact with a reservoir. This means that P e {x\, cc 2 ) 
should vanish when X\ — or L. In the dimensionless 
variables u = x/L, Eq. I)19|l becomes 



5{u x -u 2 ) (21) 



with the boundary conditions P(0,u 2 ) = P(l,u 2 ) = 
and where E^ — D/L 2 is the Thouless energy of the 
wire. The solution can be written as a sum of a special 
solution of the complete equation plus a linear combi- 
nation of the two linearly independent solutions of the 
homogeneous equation. The coefficients of this combina- 
tion can be chosen in such a way as to fulfill boundary 
conditions. The solution reads: 



P(ui,u 2 



L 



-z\u%— Tial 



2Dw 2 z 

e z(u 1 +u 2 ) _|_ e z(«i+M2-2) _ e z(«i— « 2 ) _ g— z(ux— « 2 ) 



= 2z 



1 



(22) 



where, for e > , z = ye/^tT e ~" /4 - Inserting Eq. J2U) 
into Eq. I|18[) and assuming that the size of the junction 
is small, we have 



(-PV 2 - ie)P £ (r 1 ,r 2 ) = <5 3 (n - r 2 ) 



(19) 



G^pL 



167r 3 c^unDw 2 



(23) 
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where 



Re 



sinh(zwo) sinh(z(l — uq)) 



z sinh(z) 



(24) 



and Uq indicates the position of the junction along the 
wire. Introducing the conductance of the wire G\y = 
2w 2 e 2 DvN / L and using Eq. Q15|). the expression for the 
current simplifies to 



I = 



4eG 



ir 



de 4>(e)H„ (e/2 - eV,-e/2 - eV) . (25) 



The noise obeys an identical expression with — > 
AeH+. Note that in the limit of an infinite wire the func- 
tion 4>{e) ~ diverges at low energy, as expected 
for one dimensional diffusion. Substituting Eq. <|14[) with 
Eq. I|2UI) into Eq. (|25|l one obtains current and noise for 
any temperature, voltage, and position along the wire. 
Let us consider for simplicity one specific case: zero tem- 
perature and voltage < V < U/2. It is not difficult to 
show that in this case 



2eG 



C 2 

and S = 2—^- 



7^+ ( 26 ) 



with 



4>± = (l-uo) / 4>{e)de+[{l-u Q ) 2 ±ul] / <p{e)de . (27) 

J2eV JO 

Differentiating with respect to V Eqs. (|26[) gives the 
differential Fano factor 



HeV) 



(dS/dV) 
(2edI/dV) 



(28) 



In our case T is given by the following expression 



(1 - u )<t>{2eU - 2eV) - u (2u Q - l)<t>{2eV) 
(1 - u Q )(t){2eU - 2eV) + u <j)(2eV) 



T{eV) = 2 

i i ■ a, i km in — it \ i -; //,,< n in ; 

(29) 

Note that normally the Fano factor is positive since 
it corresponds to the (positive) current noise divided by 
the absolute value of the current. The differential Fano 
factor is then usually defined in such a way that a factor 
sign(J) is implicit in the definition. This makes it usu- 
ally positive, since current and the noise are supposed 
to increase with the voltage bias. But in our case the 
wire is off-equilibrium, thus the voltage V has not the 
usual meaning of bias voltage between two systems at 
local equilibrium. Along the wire it is not even possible 
to define a chemical potential, since fermions are not at 
equilibrium. Thus the concept of difference of potential 
between the superconductor and the wire at the position 
of the contact is ill-defined. The potential V is neverthe- 
less well defined and can be used to study the evolution 
of the conductance or of the differential noise. Clearly 
the sign of the current needs not be the same as that 
of V. For this reason we use the definition l|28|l as it is, 
without changing the sign according to the direction of 



the current. This allows a more simple representation 
in Fig. |21 One should keep in mind that the change of 
sign has no special meaning in this case, since increasing 
the voltage bias V can well decrease the current and the 
noise at through the SIN junction. 

Let us discuss briefly some simple limiting cases of 
Eq. (|29|l . For uq — and 1, when the junction is at 
the extremities of the wire, Eq. (|29|) gives T = ±2. This 
is expected from Eq. Q since in this case the normal 
metal is actually at equilibrium. The sign change is sim- 
ply due to fact that for x — the potential of the wire 
U is greater than the potential of the superconducting 
point V, while for x — L the situation is reversed. For 
< i*o < 1, T is in general different from 2. An interest- 
ing point is the middle of the wire uq = 1/2. In this case 
we have: 



^(eV,uo = l/2) = 2 



(2eU - 2eV) 



4>{2eV) + (j){2eU - 2eV) 



(30) 



For V = U/2 we then obtain T = 1 exactly for any form 
of <j>(e), while for V = T should become very small 
on general grounds, since we have seen that for e — > 
the function (j>(e) is expected to diverge in the infinite 
system. 



K 




0,2 0,3 

v/u 



FIG. 2: Differential Fano factor T as a function of V/U . 
Different curves represent different values of wo going from 
to 1 in step of 0.1 starting from the top. Departure from the 
value 2 or -2 is due to the energy dependence of the tunnelling 
matrix element. 



For arbitrary values of u$ one can see a crossover be- 
tween these limiting behaviors. We report in Fig. |2Jthe 
results predicted by Eq. for U = which is a 

typical value for experiments. Wc would like to empha- 
size that the whole energy dependence seen in this plot 
is due to interference of electronic waves, since it stems 
from the e-dependence of 4>{e). 
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B. Comparison with circuit theory approach 

It can be useful to obtain current and noise with an 
alternative technique, often called circuit theory, devel- 
oped by Nazarov and coworkers The assumptions 
behind circuit theory are the same as the ones that have 
been used to arrive to Eq. (|29l) . but its scope is much 
larger, since it can be applied for arbitrary transparency 
of the interface. On the other hand its implementation is 
often numerical and can become extremely cumbersome 
in higher dimension. In the case at hand of a wire the 
problem can be solved rather easily and it is worth to 
compare the results of the two techniques. 

We follow Belzig and Nazaro^ in modelling the wire 
within semiclassical Green's function technique. The 
wire is describled by N nodes, each node represents a 
part of the wire small enough such that the inhomo- 
geneities inside can be neglected. Each node is thus 
completely described by one Usadel 4x4 matrix Q in the 
Nambu(")-Keldysh(~) space with the condition Q 2 = 1. 
The nodes are connected by tunnel junctions with con- 
ductance gx = [N + l)Gw in such a way that the total 
conductance of the wire is G\y- Kuprianov-Lukichcv 33 
boundary conditions plus the decoherence term in Us- 
adel equation give the following condition to be satisfied 
at each node i 

[A,ft]=0 (31) 

where 

A i =g T {Qi+i + Qi-i)-2i^ w n, (32) 

T3 = f 3 <S> I, and (7^, Tj(i ^=1,2,3) are Pauli matrices (a 
will appear in the following). The boundary conditions 
at the two extremities are given by the bulk solution for 
the normal metal at equilibrium: G% = f 3 g> a 3 + (/to + 
fwh) <S> (5-1 + i^a) with f L0 = 1 - /+ - /_, f T0 = /_ - 
/+, /±(e) = f(e ± eU L / R ), / is the Fermi function at 
temperature T, and Ul/r is the voltage bias at (L) x = 
and at (R) x = L. Eq. l|31|l has a simple analytical 

solution: Q = Ai/^A 2 , and the full set of equations can 
be solved by iteration leading to an explicit numerical 
value for Gi at each node. 

This describes a normal wire in good contact with two 
normal reservoirs kept at some voltage bias Ul and Ur 
with respect to (an arbitrary) ground. We add now the 
superconducting tunnel contact with conductance Gt at 
node s, with 1 < s < N. Since it is easier technically to 
keep the superconductor at zero voltage with respect to 
ground, we set Ul — U — V and Ur, — —V. The tunnel 
contact with the superconductor modifies the conditions 
(|32[) for node s in the following way: 

A s = g T (Qs+i + Qs-x) + G T Gs(x) ~ 2l J^W f * • ( 33 ) 

(Since Gt <C <7t we can neglect corrections to the Thou- 
less energy due to the presence of the superconducting 



reservoir.) The Green's function Gs(x) is the supercon- 
ducting bulk solution of the semiclassical equations, mod- 
ified by the introduction of a counting field \- This means 
that Gs(x) = e"' xfK Gs ^ xfK with G° s = f 2 ®l (since all 
energies are smaller than the superconducting gap) and 
t k = T3 (g) cti . The circuit is shown in the inset of Fig. 
Once the solution to the system of equations has been 
found it is possible to obtain the current, the noise, and, 
in principle, all higher cumulants of the current flowing 
through the SIN junction by evaluating 

J(X) = ~ J deTr{f K [G s ,Gs}} ■ (34) 

The current is given by J(x = 0) while the noise is 
proportional to its first derivative with respect to %: 
S = —2iedJ/dx- Numerically one can obtain both at 
the same time, since for very small x t ne first is propor- 
tional to the real part of J(x) while the second to its 
imaginary part. 

We have thus found numerically the current and the 
noise for U = 200E^ and different values of V. The 
numerical results presented refers to 18 nodes. The tun- 
nel conductance Gt was chosen much smaller than gT 
(gr/GT ~ 10 3 ) and we have verified that both current 
and noise scale with G T as predicted by Eq. I|26|l . 




0.1 0.2 0.3 0.4 0.5 
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FIG. 3: Fano factor as a function of V/U for x/L = 9/18, 
8/18, 4/18, and 1/18, from top to bottom. Full lines are 
obtained with Eq. 12611 . The dots indicate the result obtained 
with circuit theory. The inset shows the circuit used in the 
calculation. 

In Fig. 13 we show the Fano factor as a function of 
V/U (U = 2Q0E%, V < U/2) for different value of s, the 
node in contact with the superconductor. The agreement 
between the two approaches is rather good. 

A strong effect of the non equilibrium distribution is 
the divergence of the Fano factor for certain values of 
V or uq, while at equilibrium and at zero temperature 
Eq. predicts invariably F = 2. In the case shown in 
Fig. 0thc divergence appears for x = L/2 and V = U/2. 
This divergence can be easily understood, since at these 
values of the bias and position along the wire the current 
vanishes with = N<- ^ [cf. Eq. @] such that 
the resulting S is non vanishing [cf. Eq. i|lu|) ]. Even 
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if this is a large effect it is mainly due to the fact that 
the system is out of equilibrium, and it is thus not a 
direct signature of the importance of the interference in 
the noise. Interference appears instead directly in the 
differential Fano factor T . 

A more stringent test on the validity of both ap- 
proaches is thus to compare T{eV). Our rather simplified 
numerical approach to the circuit theory allows to extract 
the derivatives of S and / with a limited accuracy. Nev- 
ertheless the agreement of the analytical expression (|29[) 
and the numerical results shown in Fig. |I|is reasonable, 
and it remains within the numerical error of the circuit 
theory calculation. 

There are actually other reasons for possible discrep- 
ancies between the two approaches. In the analytical cal- 
culation we have not taken into account that the Fermi 
distribution changes along the wire due to the electric 
field, while this effect is included in the circuit theory 
approach. 




0.1 0.2 0.3 0.4 0.5 

V/U 



FIG. 4: Comparison of the differential Fano factor T given 
by Eq. 1291 (full lines) and by the numerical circuit theory 
simulation (dashed lines). The four dashed curves are ob- 
tained for U = 20QE% and u = 1/18, 4/18, 8/18, and 9/18 
(N = 18). 



IV. CONCLUSIONS 

In conclusion we have found a general expression that 
links the subgap noise and the current in a NIS junction. 
The interference of electron pairs leads to a non-linear 
dependence of current and noise on voltage, but their 
ratio is fixed by Eq. as long as the electrons in the 
normal metal are in equilibrium. Out of equilibrium, the 
Fano factor becomes non universal and we have computed 
it in a feasible experiment. In particular we have studied 
the differential Fano factor, which is a direct measure of 
the importance of interference of electronic waves. The 
validity of Eq. (Q) is quite general, it does not depend on 
the interactions in the normal metal, for instance. Thus 
detecting a departure from this prediction can be a strong 
experimental indication that the normal metal is out of 
equilibrium. 

From the technical point of view we have verified that 
the tunnelling approach agrees with the circuit theory 
approach. Both stem from the semiclassical theory of 
current fluctuations, but it is clear that the tunnelling 
calculation is often much simpler than solving the full 
Usadel equations. It is thus useful to verify that it can 
give reliable quantitative prediction for both the current 
and the noise in a specific example. 

Let us briefly discuss a possible experimental test of 
Eq. 0J. Due to the large barriers induced by oxides, 
experiments on shot noise in tunnel junctions are not 
as developed as they are for the transparent junctions. 
The only data available at this moment^ appear to agree 
reasonably well with our prediction. Indeed, current and 
noise both show a strong non-linear behavior, but their 
ratio follows the simple relation 

We are indebted to M. Houzet for useful discussions. 
We acknowledge financial support from CNRS/ATIP-JC 
2002 and Institut Universitaire de France. 
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